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Abstract 

The connection between the poloidal and the toroidal rotation of 
plasma in tokamak is important for the high confinement regimes, in 
particular in reactor regime. The sudden onset of closed convection 
structures in the poloidal section, due to the baroclinic production of 
vorticity, will sustain a fast increase of the poloidal velocity and a sub- 
stantial effect on the toroidal rotation. However this is limited to the 
short time of the onset transition. In real plasma however there is ran- 
dom generation and suppression of convection cells and the sequence 
of these transient events can prove able to sustain the effect on the 
toroidal rotation. We formulate a simplified model which consists of 
a laminar sheared, regular, flow situated at the boundary of a region 
of drift-wave turbulence. Vortical structures, randomly generated in 
this turbulent region are spontaneously advected toward the flow and 
are absorbed there, sustaining with their vortical content, the sheared 
flow. We examine this dynamics in the wavenumber k space, using 
reasonable approximations. We derive a system of equations (which 
is a class of Davey-Stewartson system) and find that indeed, the vor- 
tices advected and absorbed into the layer can preserve its regular, 
poloidal, flow. 



1 Introduction 

Mixed regimes consisting of large scale flows (H-mode, Internal Transport 
Barriers) and turbulence are expected to be the current state in ITER. The 
rotation, either spontaneous or induced, will play a major role in the quality 
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of the confinement. The efficiency of the sheared poloidal rotation to control 
the instabilities is considerably higher than that of the toroidal rotation but 
it is usually assumed that the poloidal rotation should be at the neoclassi- 
cal level due to the damping induced by magnetic pumping. This is true 
if the drive of the poloidal rotation (besides the neoclassical one) relies on 
the Reynolds stress produced by a poloidally cuasi-symmetric turbulence. 
Much higher drive of the poloidal rotation is however provided by flows as- 
sociated with the convective structures that can be generated in the plasma 
cross- section beyond a threshold in the plasma pressure gradient. This drive 
overcomes the damping due to the magnetic pumping and the poloidal ro- 
tation is sustained. Cells of convection consisting of closed, large scale flows 
with unique sense of rotation can be spontaneously generated, triggered by 
streamers sustained by the baroclinic term able to generate vorticity. Similar 
to the Reyleigh-Benard first bifurcation (from purely conductive to convec- 
tive regime), the onset is very fast and the drive exerted on the poloidal rota- 
tion leads to a fast time variation of the polarization radial electric field. This 
is sufficient to create a distinction between the phases (first half, second half) 
of bounce on a banana of trapped ions and, implicitly, leads to acceleration 
in the toroidal direction. This phenomenological model has been proposed 
recently [1] and the estimations are compatible with the observed effect of 
reversal of the toroidal rotation in tokamak. There are several aspects that 
still must be investigated and in particular the nature of sudden-impulse of 
the generation of the poloidal flow as envelope of the peripheral velocity of 
the convective structures. It is only on a very short time interval (onset of 
convection) that the time variation of the radial electric field can produce a 
toroidal effect on the bananas. If the convection cells become saturated and 
stationary then the efficiency of sustaining toroidal acceleration of bananas 
decreases rapidly to zero. The process must be repeated if it is to sustain its 
effect on the toroidal rotation. 

On the other hand, it is highly improbable that the structure consisting 
of cells of convection breaking the azimuthal symmetry can exist as a station- 
ary state. This is because the convection in the cells is much more efficient 
as radial transport of energy than any diffusive process sustained by small 
scale turbulence and the equilibrium profiles would be modified, suppressing 
the origin of convection. It is more realistic to regard the generation of con- 
vective structures as a stochastic process consisting of a random sequence of 
transient events, taking place on a range of spatial scales between p s -scale 
drift vortices up to cells of convection with dimensions comparable with the 
minor radius a. We should study the way these transitory structures create 
the impulsive increase of the poloidal velocity such as, on the average, to 
sustain the toroidal effect on bananas. A marginal stability regime can be 
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established, where the external sources sustain the gradient of the pressure, 
the convection rolls are formed transitory and the excess of the gradients rel- 
ative to a threshold is transformed into pulses of heat and momentum in the 
radial direction. This problem requires a treatment based on the stochastic 
nature of the generation and destruction of the transversal convection cells, 
seen as a random transitory process. Such a treatment simply does not exist, 
not even for the famous Rayleigh-Benard system. 

This problem is important per se, beyond the phenomenological model for 
the reversal of the toroidal rotation. This process has been identified in exper- 
iments on a linear machine [2] and explained in the context of residual stress 
able to sustain the coherent flow. The sequence consisting of tilt, stretch 
and absorption of drift turbulent eddies was found to contribute significantly 
to the Reynolds stress. A similar physical picture has been proposed by 
Spineanu, Vlad, Itoh and Itoh JBJ to explain the formation of sheared flow 
layers (acting as transport barriers) via tilting instability, in analogy with 
the winds of the Rayleigh-Benard second bifurcation. The absorption of a 
coherent convection cell inside a layer of sheared rotation is mentioned by 
Drake et al. jl]. It has also been discussed in similar terms in Refs.JS], [5J. 

The process in which strongly tilted cells of a convection pattern coalesce 
to form a contiguous layer of poloidal flow is the final phase of the instability 
described by Rosenbluth and Shapiro [7] . The final state consisting of layers 
with alternating direction of flow is favoured by the fact that this pattern is 
the exact solution of an equation representing the extremum of an energy- 
like functional [8]. However, in order to incorporate the effect of vortical 
structures having a wide range of spatial scales that are absorbed into a pre- 
existing layer of sheared flow one needs a different approach, not elaborated 
to date. It may be useful to compare with the generation of topological 
(global) vortices in a rotating Bose - Einstein condensate [9]. When the 
angular momentum of the condensate exceeds a threshold a new vortex is 
generated at the edge and moves toward the center. Our problem is the 
opposite process, in which vortices or convection structures accede to the 
layer of coherent flow and sustain it with their vorticity content. 

The structure of the paper is as follows. We formulate a model which 
is relevant for the physical problem of sustainement of poloidal rotation by 
stochastic rise and decay of vortices. We estimate the average density of high 
amplitude vortices that are generated inside a volume of plasma (in our case 
it is question of an area within the poloidal cross-section since the third - 
toroidal - axis is considered irrelevant). Then we remind the known mecha- 
nisms on motion of a vortex on a background of gradient of vorticity, in order 
to see how random elements of vorticity will join the flow layer where they 
will absorbed. Finally, we examine this dynamics in the wavenumber space 
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k and derive a model system of equations reflecting the interaction between 
the laminar flow and the incoming drops of vorticity. This is actually the 
Davey Stewartson system and indeed, the class of solutions called "solitoffs" 
show that the final state exhibits unidirectional (poloidal) alignment of the 
resulting flow. This confirms that the elements of vorticity coming from the 
turbulent region into the laminar flow do not perturb but are conforming to 
the geometry of the flow, i.e. in the poloidal direction. 

2 Formulation of the model 

In connection with our problem of stochastic stationarity of sustainment of 
poloidal flow by intermittent rise and decay of convection cells we construct 
a simple model to be investigated. We consider a layer of laminar sheared 
rotation where we assume that the instabilities are suppressed (this is an 
ideal representation of the zonal flow or of the if- mode layer). The layer is 
connex to a region where drift-type turbulence exists. In the drift-turbulent 
region there is stochastic generation of strong vortical structures, since these 
are a possible solution of the drift type instability. These are not solitons but 
are very robust, cuasi- coherent structures characterized by a high magnitude 
of the vorticity. 

Since the background has a gradient of vorticity a vortex will move in a 
direction which depends on the sign of its vorticity relative to the one of the 
background: the clumps (positive sign circulation vortices) are ascending 
the gradient of vorticity, toward the maximum of vorticity in the rotation 
layer and the holes (negative sign circulation vortices) are moving toward 
the minimum of vorticity, equivalently - they are repelled from the sheared 
layer [TO]. In short, vortices of only one sign are joining the layer. The 
clump vortices will be absorbed by the sheared velocity layer and they will 
contribute with their vorticity content to the momentum of the background 
flow. This is a source of momentum which sustains the sheared flow, in 
particular sustains the sheared poloidal rotation against the damping due to 
the magnetic pumping. 

We take the layer of plasma flow in the y (= longitudinal, streamwise, 
poloidal in tokamak) direction with sheared velocity varying in the x (= 
transversal on the layer of flow, radial in tokamak) direction. It is described 
by a scalar function A (x, y, t) which is the complex amplitude of slow spatial 
variation of the envelope of the eddies of drift turbulence. From A (x, y, t) the 
components of the velocity are derived as v s = — VA (x,y) x n/S. We also 
define a scalar function <p (x, y, t) which represents the field associated to the 
incoming vortices that will interact and will be absorbed into the sheared, 
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laminar, flow, transferring to it their vorticity content. The process in which 
a localized vortex is peeled-off and absorbed into the sheared flow may be 
seen as a transfer of energy in the spectrum from the small spatial scales 
involved in the structure of the vortex toward larger spatial scales of the 
sheared flow, with much weaker variations on the poloidal direction. This 
can be seen as a propagation in k space. 

We must first estimate what is the number of high amplitude, robust vor- 
tices that are generated by the drift-turbulence in a unit area. This problem 
is treated in the Appendix 1 where it is shown that the density of strong vor- 
tices is n = 3/ (vrA 2 ) where A is the typical size of the drift eddies normalized 
to p s . 

3 How these vortices ascend to the flow layer? 

The answer is: due to the existence of the background gradient of vorticity. 
This has been clearly described for non-neutral plasma in [TU] and further 
applied in the physics of the atmosphere and of protoplanetary disks. 

Since the background has a gradient of vorticity a vortex will move in a 
direction which depends on the sign of its vorticity relative to the one of the 
background: 

• the prograde vortices are moving toward the maximum of vorticity and 

• the retrograde vortices are moving toward the minimum of vorticity. 

The prograde vortices eventually are absorbed by the sheared layer and 
they contribute with their vorticity content to the momentum of the back- 
ground flow. This is a source of momentum which sustains the sheared flow, 
in particular sustains the sheared poloidal rotation against the damping due 
to the magnetic pumping. 

A clump of vorticity (positive circulation T v > 0) goes up the gradient 
of background vorticity. For a sheared velocity flow with poloidal velocity vg 
increasing linearly with radius, with a coefficient Q 

v g = Qr 

but in a layer which has background vorticity (' Q along r , the radial motion 
of the positive vortex (position r v ) is 

^=Q 2 ln (^)arctan(|Q|t) 
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where I = y/\T v / (2ttQ)\ [II]. 

The negative-circulation vortices, < 0) called holes are moving in the 
opposite direction therefore they remain in the turbulent region and are re- 
pelled by the layer of sheared flow. This means that only vortices of a certain 
sign ascend to the flow and are absorbed in it. In the following we will anal- 
yse this process as an interaction between the two fields: A (x, y, t) and the 
field representing the incoming vortices. 

4 Dynamics in the k - space of the interaction 
between the laminar flow field and the field 
of the vortical cells 

In order to derive a dynamic equation for the scalar streamfunction of the 
sheared flow A (x, y, t) we remind the properties resulting from multiple 
space-time analysis of a turbulent field, like the drift waves interacting with 
a coherent, sheared flow. The typical equation in the hydrodynamic de- 
scription is Rayleigh-Kuo or barotropic equation and the multiple space-time 
scale analysis leads to the Nonlinear Schrodinger Equation for the amplitude 
W(x,y,t) of the envelope of the oscillatory waves excited in the fluid |12j . 
The general structure of this equation is 

i ?^^ = AW(x,y,t)-j\W\ 2 W* (1) 

The interesting term is the nonlinearity (the last term) which may be seen 
as arising from a self-interaction potential V [W (x, y, t)] of (p A -type. The 
later potential results naturally from the multiple scale analysis and it is 
the effect on the envelope amplitude of the nonlinearity consisting of con- 
vection of the vorticity by its own velocity field. We must work with an 
amplitude A (k x , k y , t) in the spectral k— space, which we see as the envelope 
of oscillations representing propagation of perturbation in the spectrum, this 
propagation being a reflection of what is the process of the incoming vortex 
dissolving itself and being absorbed by the flow, - in the physical space. A 
vortex of small spatial extension is a concentration of energy at high k in the 
spectrum and its dissolution into the flow means that this energy leaves the 
high k part to go to smaller k part in the spectrum, since small k is specific to 
poloidal coherent flow. This is a propagation of energy in k space. Looking 
for an analytical model for A in the k-space we will use the structure which is 
suggested above. Although there is no rigorous proof for this, it is suggested 
by the similar processes related to the Landau-Ginsburg nonlinearity [9]. 
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We expect to have a diffusion in k— space and the same self-limiting non- 
linear term 



. dA(k x ,ky,t) _ (<P__&\ 



m \dk 2 x + dki) 



A (k x , k y , t) - 7 \ A\ 2 A* + S 



(2) 



and we have added a source that represents the interaction of the basic flow 
A with the vortices that are randomly generated in the drift turbulence and 
are joining the flow. 

The addition of "drops" of vorticity is represented by the k-space scalar 
function 



Several space scales are present in the structure of the function (x, y, t) in 
real space. The reason is that we must consider that the vortices that are 
continuously generated in the turbulent region have two characteristics: 

• the space extension of individual vortices coming into the layer has a 
wide range of values. The range extends from vortices of the dimension 
p s (the natural result of excitations sustained by the ion-polarization 
drift nonlinearity) up to convective events representing streamers with 
cuasi-closed (roll-type) geometry and spatial extension given by the 
inverse of kg = m/r for m of only few units. 

• the spatial distribution of the place where these vortical structures are 
reaching the layer of sheared flow is random over practically all the 
circumference of the poloidal flow. 

We will draw conclusions about the k-space profile of based on the 
contributions coming in from various scales in the real space. The most sim- 
ple description neglects the random aspect of the content of and simply 
enumerates the possible spatial scales. For this we adopt a discrete represen- 
tation of the dimensions of vortices, starting from a minimal size of the order 
of p s up to a low-m fraction of the circumference 2nr of the poloidal rotation 
layer, this last scale corresponding to large convection roll-type events. 

The above elementary assumption, i.e. the discrete number of space ex- 
tension of vortices, implies that the function (k x , k y ) presents higher ampli- 
tudes around wavenumbers k which correspond to the inverses of the spatial 
scales. The simplifying assumption which takes uniform spatial distribution 
of the events of arrival of vortices (of any scale) to the layer implies period- 
icities on various spatial scales. 

We consider separately the y direction (poloidal). One periodicity can be 
at the level of a smallest vortex, Jty\ This means that a sequence of maxima 



{k x , ky, t) 



(3) 
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of (x, y) appear with spatial periodicity 2ix/k^ and this corresponds to the 
addition of smallest vortices into the flow layer. Further, on a longer space 
scale there is another periodicity, k v 2 \ This means that the sequence on 
the lower scale (with much denser spatial granulations of vorticity, 2-ir/ky 1 ^) 

(2) 

is modulated on a longer space scale with a periodicity 2i\jk y . This is 
because larger vortices (but not yet convection rolls) are absorbed into the 
flow layer. Large vortices carry smaller vortices. These considerations can 
be extended to other, intermediate, levels of periodicity. Finally we can have 
the largest (accessible) spatial scale 27r/k y n ^ where a convective event occurs 
as a result of a stochastic event, or a streamer sustained by the baroclinic 
term and having the close-up property. 

We note that the periodicity in the poloidal direction y cannot occur 
without a corresponding periodicity in the x = radial direction, as far as the 
vortices and convective events are concerned. It is reasonable to assume a 
linear dependence : small scale periodicities in the poloidal direction (small 
scale robust drift vortices) occur on a similar radial scale, i.e. again small 
periodicity in the x = radial direction. With larger y wavelengths, we expect 
to involve also larger radial wavelengths since the physical process involves 
more fluid motion in the closed-up convection cell. 

This suggests a proportionality between the periodicities on the poloidal 
direction, y, and the periodicities on the radial, x, direction, in the struc- 
ture of (f>(k x ,ky). The fact that there is (assumed) proportionality on the 
periodicities on y and on x, requires for analytical description a hyperbolic 
operator 

/ d 2 1 d 2 \ 

where the D 'Alamebertian operator Ok x ky = Jj3 — ^ * s the signature of 
the correlated periodicities on the poloidal and radial directions of the field 
representing the incoming drift vortices - up to convection rolls, generated 
in the drift turbulence region, moving against the gradient of vorticity, and 
finally absorbed into the rotation layer. The "velocity " a is the ratio of the 
periods on the two /c-space directions and we take it a ~ 1. The centers of 
the vortices reach points that belong to the lines y = ±ax. The general form 
of the solution of the homogeneous equation is 

(f)(k x ,k y ) ~ exp(ik x x + ik y y) (5) 

When is "free" it is double periodic and it verifies the above equation. 
However the destruction (peeling-off) of the incoming vortices and their ab- 
sorption through interaction with the background flow, A, is the loss of the 
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double periodicity. The departure relative to the k-space double periodicity 
(the "free" state of <fi) feeds the flow field A (k x , k y ,t) such as to enhance it 
in spectral regions comparable to those where the vortices are localized, i. e. 
higher |k|. Basically, if the free A (k x , k y , t) is localized on k y close to k y ~ 
(equivalent to an uniform poloidal flow), the source induced by is a charge 
which acts through a Poisson equation on the "potential" given by the am- 
plitude AA*. We should recall that we want to represent the physical process 
in which incoming vortices transfer via absorption their vortical content to 
the sheared flow, sustaining and/or enhancing the poloidal, y, rotation. This 
is mainly manifested as the evolution of an initial oscillation generated by </> 
on the fcy-spectrum with collection of the energy in a region close to k y ~ 0, 
since this corresponds in real space to the uniformization of the flow in the 
poloidal direction, after an event of absorption of a vortex. Then, ignor- 
ing the less significant periodicity in k x direction, corresponding to periodic 
perturbation propagating in the x (radial) direction, the interaction may be 
represented as 

The constant /3 is a measure of the permittivity of the equivalent electro- 
static problem. Using the specific terminology measures the decrease of 
the effectiveness of the effect of the charge (right hand side of the Eq.([n])) on 
the potential |v4| 2 due to the "polarization". Essentially, the amount of de- 
struction of the incoming vortex, measured by the departure from the "free" 
state by the right hand side of the Eq.([6]) is transferred to the background 
flow by first exciting elementary waves in real space, while A is only an en- 
velope. The perturbation of the sheared flow A (via the left hand side of the 
Eq.([n])) results from the nonlinear interaction of these excited waves, i.e. by 
convecting the modified vorticity with the background velocity flow. This is 
the last term of the left-hand side of Eq.([2]). Then it is reasonable to assume 
w 7. 

Returning to the equation for A (k x , k y ,t) we specify the interaction be- 
tween the two fields A and <ft m the simplest way 

S = 2(/)(k x ,k y ,t)A(k x ,k y ,t) (7) 

Then 

dA(k x ,k y ,t) ( d 2 d 2 \ A „ , , . A . „ 

' at = (*I + eiq) A fc »" () - 7 W A + 24,A (8) 

to which we add the equation for (ft, 

d 2 d 2 \ d 2 

0(/^,t)«- 7 _|A| 2 (9) 



dkl dkl ) T v ' y ' ' 'dk 2 
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These two equations are known as Davey-Stewartson system and our 
case is DS-I [T3]. It is exactly integrable and several analytical solutions are 
available, in terms of Riemann theta functions or Jacobi elliptic functions. 
Of particular importance for our physical problem is the long wave limit 
solutions, as obtained by Chow using the Hirota method. We reproduce here 
his result [T3] 

. r(l-k) 



tanh (sk x ) — y/ksn (rk y , k) 
1 + v^sn (rky, k) tanh (sk x ) 



exp (-iQt) (10) 



2 k (ksn 2 (rky, k) + tanh 2 (sk x )) + (1 + k 2 ) -y/fcsn (rk y , k) tanh (sk x ) 
(l + v^sn (rky, k) tanh (sk x )j 

where s = r (1 + k) /2, ft = 4r 2 (1 - E/K)-r 2 (1 - kf /2, K and E are the 
complete elliptic integrals. These solutions are called solitoffs and they decay 
in all directions except a preferred one. They are " semi- infinite solitary 
waves". For example, for k — 0.7 and r = 1 the solution shows a peak 
of \A\ around k y » which in physical terms means quasi-uniform flow 
along the poloidal direction. This is compatible with what we expect from 
this analytical model: the turbulence advects "drops" of vorticity into the 
rotation layer and this does not destroy the flow but sustains it. It is actually 
a discrete Reynolds stress, or a process which is reversed version of the usual 
Kelvin-Helmholtz instability. 



5 Conclusion 

In conclusion, this simple analytical analysis supports the idea that the ran- 
domly generated vortical motions, including the convection rolls that rise and 
decay transiently in the plane transversal to the magnetic field, are able to 
sustain poloidal rotation. When the time scale is very short, as for convec- 
tion rolls driven by the baroclinic effect, the poloidal flow arising as envelope 
of the rolls has fast time variation and the neoclassical polarization effect 
acts effectively on the toroidal rotation. This process accompanies every 
event consisting of configuration of a closed convection pattern. Then, even 
if these events are random, their statistical average is an effective source of 
poloidal flow. 
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A Appendix. Estimation of the number of 
vortices of high amplitude that are gen- 
erated in every unit area of a region of 
drift-wave turbulence in plasma 

In this Appendix we discuss the problem of estimating the average number 
of strong vortices that are generated in a two-dimensional region of drift- 
turbulent plasma. Being a complex problem, it will be the subject of a 
separate work. We limit ourselves here to mention the main steps of this 
calculation, revealing the necessity to link apparently different systems. 

We adopt the definition that a strong, robust vortex of high amplitude is 
one which is a physical realisation of the mathematical singular vortex, the 
latter being characterised by u> (r) ~ 5 (r). This approximation is useful since 
we will now look for the statistically averaged density of singular vortices in 
a region of turbulence. 

The following connections can be made: 

1. The drift wave turbulence in the hydrodynamic regime and at spatial 
scales that are much larger than p s is governed by the cuasi-three- 
dimensional Hasegawa-Mima equation. However, for strong magnetic 
field, there is not too much difference to ignore the intrinsic length, the 
Larmor radius, and to consider the approximate description given by 
the scale-free hydrodynamic model, i.e. the 2D Euler equation 

d 

—Aip + [(- Vtp x e z ) • V] Aip = 

2. It is known that this equation is equivalent to a model of point-like 
vortices interacting in plane by a self-generated long-range (Coulom- 
bian) potential as known from prestigious theoretical work Kirchhoff, 
Onsager, Montgomery, etc. (see [13]). 

3. The dynamics of the point-like vortices has been mapped onto a field 
theoretical model consisting of a complex scalar (matter) field </>, <j>\ 
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and a gauge field (which mediates the interaction) [T5]. The La- 
grangian density consists of the kinetic term for <f), the Chern-Simons 
term for N 1 and a nonlinear interaction. The field variables are in the 
algebra su (2) due to the vortical nature of the physical elements 

with E± the ladder generators in the cartan algebra. Introducing p x = 
|0 1 | 2 , p 2 = |0 2 1 2 , the vorticity is 

w = — (px - p 2 ) 

and in the relaxed states (self-dual) we have 

4 

Aip H — sinh ip = 

«; 

here ■?/> is the streamfunction of the flow, with ip — In p, p t — p^ 1 = p. 

4. We have noted that this equation is exactly the equation describing the 
surfaces in real 3D space with the property of constant mean curvature 
(CMC). If the curvatures of the surface are Ki and k,2 the CMC surfaces 
have 

H = Ki + K 2 = 1 

5. The fact that the self-duality states and the CMC surfaces are gov- 
erned by the same equation sm/i-Poisson has a deep significance. Ex- 
tending the parallel between these two system at states which are not 
at self-duality or at CMC, {i.e. the dynamical states preceeding the 
asymptotic relaxed states) we have established the following mapping 

(«i + K 2 f 
Pi = 

Ki - K 2 
p 2 = Ki - K 2 

6. It then results that the umbilic points of a surface in real 3D space, 
i.e. a point where 

Ki = K 2 

is a point where 

P X — > CO 

P2 ->• 

and this means that 

\oj\ —> oo 

at these points. 
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7. Now we see that looking for the points where there is a very high value 
of the vorticity (almost singular vortices \u\ — > oo) means to look for 
the umbilic points (k± = K2) of a surface. 

8. We actually consider a turbulent region where the strong, cuasi- singular 
vortices appear at random. Equivalently, we must consider a surface 
in 3D which has fluctuating shape. The statistical properties of the 
random strong vortical structures arising in turbulence is the same as 
the statistical properties of the set of umbilic points on the fluctuating 
surface. This statistics is taken to be Gaussian. 

9. We then have to find the statistical properties of the set of umbilic 
points on a surface whose shape fluctuates with the Gaussian statistics. 

This problem is solved in Ref. [16]. For a Gaussian correlation of the 
surface heigth p ~ exp (— r 2 /A 2 ), the density of umbilic points is 

3 

n = — o 

and this is also the density of high amplitude robust drift-type, vortices 
generated by turbulence. Here A is the linear dimension of a typical eddy of 
the drift wave turbulence, normalized to p s . 
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